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o : 

' Abstract. Using the lag-luminosity relation and various BATSE catalogs we create a large catalog 

O \ of burst redshifts, peak luminosities and emitted energies. These catalogs permit us to evaluate the 

D . lag-luminosity relation, and to study the burst energy distribution. We find that this distribution can 

' be described as a power law with an index of a = 1 .76 ± 0.05 (95% confidence), close to the a = 2 

\ predicted by the original quasi-universal jet model. 

^ : INTRODUCTION 

^ : 

^ ■ Jet models predict the distribution of the isotropic-equivalent energy quasi- 

^ ! universal jet profile models predict an approximate power law distribution with index 

a = 2, where N{Eiso) °^ Ej^^[l]. The isotropic-equivalent energy Ejso is the total energy 
radiated if the observed flux were radiated isotropically. To study the distribution of 
burst intensities, we used the lag-luminosity relationship to create a burst database with 
redshifts, peak luminosities, and burst energies, and then we fit energy distributions to 
the burst database. Of course, this database can be used for other studies. 

In the lag-luminosity relation[2] the peak bolometric luminosity Lg is a function of the 
^ ', lag Tfi between two energy bands in the burst's frame — Lg = But Tq is measured 

I in our frame. We model Tg = (1 +zYto: time dilation contributes -1 to c, while the 

^ ■ redshifting of temporal structure with a smaller lag from higher energy contributes ~l/3 

^ ! (pulses are narrower at high energy). The peak bolometric luminosity is related to the 

H I peak bolometric energy flux Fb = Lb/[AkD]], where Dl is the luminosity distance. The 

peak bolometric energy flux is related to the peak photon flux P (integrated over an 
energy band, e.g., 50-300 keV for BATSE data): Fg = {E)P. The result is an implicit 
equation that must be solved for each burst: 

P = Q{{\+zrxQ)/[{E)AnD]\ (1) 

After solving eq. 1 for the redshift, Lg and Eiso can be calculated from Fb and the energy 
fluence, respectively. 

The original lag-luminosity relation was a single power law, e.g., Lg 'ig But 
this power law over-predicts the luminosity of GRB980425 (assuming this burst was 
SN1998bw). Consequently Salmonson[3] and Norris[4] suggested breaking the single 
power law; for Tg >0.35 s the power law index is -4.7. A population of nearby, long lag 
bursts resulted. 
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With a database of bursts with Eiso we can now calculate the energy distribution. 
The methodology presented here[5] can also be applied to the luminosity function. 
The probability of detecting a given energy is truncated by the detection threshold: 
p{Eiso\Eiso,ihM{a)) where Eiso,th is the threshold value of Eiso for that burst and M{d) 
is the model (e.g., the functional form of the energy distribution) with parameters a. For 
the ensemble of bursts the probability of detecting bursts with the observed energies is 

A = fjp {Eisci I Ei,o,th,iMi^)) (2) 

i 

where the product is over each burst. This probability is the likelihood for the model 
M{a). In frequentist statistics, we maximize A with respect to the parameters a to get a 
best fit value. In Bayesian statistics the likelihood is a factor in the "posterior," which 
can be used for confidence ranges and best fit values; the Bayesian approach allows the 
use of "priors" reflecting our expectations for a. 

Note that to study the energy distribution we do NOT need a complete sample in 
terms of observed fluences, only a sample that has no bias on the intrinsic Eiso. There 
can be gaps in the distribution of peak fluxes, but Eiso has to be drawn uniformly from 
p{Eiso\Eiso,thM{a)) in our sample. On the other hand, if we want the burst rate per 
comoving volume as a function of redshift, then we do need a complete sample. 

But is the resulting energy distribution a good representation of the data? The like- 
lihood (frequentist approach) or posterior (Bayesian approach) can be used to compare 
models (functional forms), but do not indicate "goodness-of-fit." However, our method- 
ology assumes the energies are drawn uniformly from p{Eiso \Eiso,thM{a)). The cumu- 
lative distribution of p{Eiso \ Eiso^th^i^)) should therefore be a straight line, and the av- 
erage value should be 1/2, with a statistical uncertainty of [12N]~^/'^ for N bursts. 

RESULTS 

We started with 1438 BATSE bursts for which we calculated lags. Of these, 1218 
have positive lags. These bursts also have hardness ratios, peak fluxes and durations. 
To calculate the average energy (E) we used the "GRB" spectral fits of Mallozzi et 
al.[6] to the peaks of 580 of these bursts. For the 858 bursts without fits we assumed 
average spectral indices a = -0.8 and j8 = -2.3. Plotting HR32 (the 100-300 keV to 
50-100 keV hardness ratio) vs. Ep shows a clear correlation which can be approximated 
by Ep=24Q HR32 keV; we used this relation for the bursts without spectral fits. 

Redshifts were calculated for this database for both the original simple power law lag- 
luminosity relation and the broken power law Salmonson[3] and Norris[4] introduced 
to incorporate GRB980425. As expected, the difference in the lag-luminosity relations 
is apparent at low redshifts: the broken power law results in a population of nearby 
bursts. There were few physically implausible high z bursts (e.g., z > 20) and thus no 
additional cutoffs on the lag-luminosity relation are required. In the absence of additional 
information, the choice between the two lag-luminosity relations depends on whether 
GRB980425 is considered to be a typical low luminosity burst. For the remainder of this 
analysis we use a single power law lag-luminosity relation. 
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FIGURE 1. Scatter plot of the isotropic equivalent energy Eiso vs. the detection threshold. 



As an aside, we found that the redshift calculation is sensitive to the value of {E). 
Calculating this quantity inconsistently can introduce errors into the resulting database. 

We calculated the energy Eiso for each burst from the redshift and energy fluence. The 
results were reasonable (see Fig. 1): few bursts had Eis„ > 10^^ erg or Ei^o < 10^^ erg. 

The energy detection threshold Eiso^th can be calculated by scaling Ei^o by the ratio of 
the threshold peak photon flux to the observed peak flux. BATSE's threshold peak flux 
was Pmin ~0.3 ph cm~^ s""^; however, the number of bursts in our sample with Eiso just 
above the threshold is suspiciously low (see Fig. 1), suggesting that the sample's true 
threshold was greater than 0.3 ph cm~^ s~^ Consequently we used Pmin ~0.5 ph cm~^ 
s~^ as the threshold, deleting bursts with P < 0.5 ph cm~^ s~^ 

The left hand side of Fig. 2 shows the likelihood surface for our sample assuming a 
power law functional form, where the two parameters are the low energy cutoff E2 and 
the power law index a (i.e., N{Eiso) °^ E^^^ for Eis„ > E2). The likelihood is maximized 
by E2 equal to the lowest observed value Eiso, although lower values are not ruled out. 
The best fit spectral index is a = 1.76 ±0.05 (95% confidence). Although (P(> Eiso))= 
0.4642±0.0089 (A^=1054, assuming P^;„ ~0.5 ph cm'^ s'^) deviates from 1/2 by 4a, 
considering the possible systematic errors (e.g., in the estimation of Ep from the hardness 
ratio), this value of {P{> Eiso)) indicates that a power law energy distribution is a fairly 
good characterization of the data. 

We also tried a lognormal energy distribution (right hand side of Fig. 2). The max- 
imum likelihood occurs at £'iso,cen = 3 x 10^^ ergs and cj£ = 2.7. The surface's shape 
indicates that the data permit a high central value of Eiso and a narrow distribution, or 
a low central value of E^so and a broad distribution. The observational cutoff truncates 
the true energy distribution, and the low energy extent is relatively unknown. We find 
(P(> £))=0.4821±0.0089 (A^=1054, assuming ~0.5 ph cm'^ s'^), consistent with 
1/2 at the 2a level. 
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FIGURE 2. Contour plots of the likelihood surface for a power law energy distribution (left) and 
lognormal energy distribution (right). The power law has a low energy cutoff £2 and power law index 
a; the contours are spaced by Alog(likelihood)=l. The lognormal distribution has a central value ^isojcen 
and a logarithmic width Oe\ the contours are spaced by Alog(UkeUhood)=10. 



IMPLICATIONS 

A quasi-universal jet profile that is a power law in the off-axis angle Q — the energy per 
solid angle e(0) oc — results in a power law energy distribution (or luminosity func- 
tion) with index a = X—ljk (hence a = 2 for k = —2), while a Gaussian profile results 
in a = 1. Lloyd-Ronning et al.[l] found that if the profile parameters are distributions, 
the luminosity functions could be approximated by power laws with ct ~ 2 for power 
law profiles and a ~ 1 for Gaussian profiles, but with curvature. The additional degrees 
of freedom introduced by varying the parameters give the jet models the freedom to fit 
a wide variety of energy distribution shapes. We find that our burst data can be fit by 
a power law energy distribution with a — 1.76±0.05 (95% confidence); considering 
only the statistical uncertainty the power law distribution is formally not a good fit, but 
with the likely systematic uncertainties the power law distribution is probably a good de- 
scription of the data. While our power law fit is inconsistent with the original jet profile 
model (k — —2 and therefore a — 2), it is consistent with the jet profile models where 
parameters are permitted to vary. 

A log-normal energy distribution also describes the data; the data permit a smaller 
average energy if the distribution is wider. 
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